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ABSTRACT 


Seven  terms  of  the  low  temperature  partition 
function  for  the  face-centered  cubic  antiferromagnetic 
lattice  with  Ising  interactions  have  been  calculated 
by  counting  bi-coloured  graphical  configurations. 

Low  temperature  series  expansions  for  the  specific 
heat.,  internal  energy,,  and  the  zero  field  magnetic 
susceptibility  are  derived.  The  transition  temperature 
kTc/J  is  estimated  to  be  about  1.9  *  .1*  obtained  by 
extrapolating  various  series  by  Pade  approximant  and 
ratio  methods.  The  specific  heat  curve  was  fitted 
to  a  logarithmic  singularity  at  T  and  from  the 
resulting  curve  an  estimate  of  the  critical  entropy 
was  obtained  by  direct  integration  of  the  area  under  the 
curve  .We  obtain  Sc  =  .312  k. 


Introduction 


1.1  Cooperative  Phenomena  in  Crystals 

Macroscopic  assemblies  can  be  divided  into  two 
categories.  In  the  first  category ,  the  microscopic  subsystems 
of  the  assembly  can  be  regarded  as  independent  and  non  - 
interacting.  Here  the  thermodynamic  properties  of  the  total 
system  are  determined  by  the  energy  levels  of  the  micro¬ 
scopic  subsystems  and  usually  exhibit  smooth  and  continuous 
behaviour.  The  second  category  includes  systems  in  which 
the  microscopic  subsystems  interact  strongly.  Here  we  find 
the  phenomenon  of  phase  transitions  occuring,  such  as 
gas  condensation,  order-disorder  transitions  in  alloys  and 
ferromagnetism . 

In  1925  Ising^1^  proposed  a  simple  model  of  such 
a  system  in  which 

(1)  All  atoms  were  identical  but  to  each  one  there 
corresponds  only  two  possible  states. 

(2)  Only  nearest  neighbour  (n.n.)  interactions 
are  significant. 

He  associated  with  each  lattice  point  a  spin 
coordinate  a  which  could  take  on  one  of  two  possible  values, 
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+1  for  spin  up  and  -1  for  spin  down.  The  interaction 
energy  between  the  j  and  kth  lattice  points  was  pos¬ 
tulated  to  be 


-Jo  .a 


j  k 


if  j  and  k  are  n.n. (nearest 

neighbour ) 


( 0  otherwise 

Thus  the  energy  is  -J  if  the  n.n.  atoms  have  parallel 
spins  and  +J  if  they  have  antiparallel  spins. 

In  addition,  if  there  is  an  external  magnetic 
field,  H,  present  each  spin  is  coupled  to  the  field  by 
and  interaction  energy  -m  H  if  the  spin  is  aligned  with 
the  field  and  +m  H  if  the  spin  opposes  the  field, 
m  is  the  magnetic  moment  per  spin. 

We  have  then  for  the  total  energy  in  a  lattice 


of  N  sites 

N  N 

E  =  -J  y  a. a.  -  H  m  )  a.  (1) 

i , J = 1  1  J  i=l  1 

j<i 

and  the  partition  function  is 


a  i  =»±  1 


L-I 

2='-  1  °N 


exp 


J 

"  ’• 

.  kT 


l 


R. 


i  J=1 


from  which  in  principle  all  the  thermodynamic  properties 


can  be  derived. 
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3. 


1.2  Antiferromagnetism  and  the  Face-centered  Cubic  Lattice 

When  J  Is  negative  we  say  that  the  system  is 
antiferromagnetic . 

In  zero  magnetic  field  the  lowest  energy  states 
for  a  ferromagnet  will  be  the  two  states  in  which  all  the 
spins  are  aligned.  However,  in  the  case  of  an  anti¬ 
ferromagnetic  lattice  the  lowest  energy  states  will  be  the  two 
in  which  the  n.n.  spins  tend  to  be  antiparallel. 

For  some  lattices  such  as  the  simple  quadratic 
(sq),  simple  cubic  (sc)  and  body  centered  cubic  (bcc) 
perfect  antiferromagnetic  ordering  is  possible  (all  spin 
up  atoms  are  n.n.  of  spin  down  atoms  only  and  vice  versa). 
Clearly  these  "loose  packed"  lattices  can  be  subdivided 
into  2  sublattices;  one  containing  only  "spin  up"  atoms 
and  one  only  "spin  down"  atoms  and  nearest  neighbour 
bonds  exist  between  atoms  in  different  sublattices  only. 

For  other  lattices  including  the  triangular  and 
f.c.c.  such  A.F.  ordering  is  not  possible  since  there  exist 
closed  paths  in  the  lattice  with  an  odd  number  of  nearest 
neighbour  bonds.  Very  little  information  is  known  about 
the  thermodynamic  properties  of  these  lattices,  especially 
in  the  low  temperature  region. 

This  thesis  is  concerned  with  the  properties  of 
the  spin  1/2  Ising  model  of  thef.c.c..  ant  if erromagnet . 
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4. 


Among  the  various  lattices  that  have  been 

studied  extensively,  the  f.c.c.  antif erromagnet  seems 

to  be  unique  in  that  it  is  the  only  close  packed  lattice 

that  has  been  shown  to  have  an  order-disorder  transition 

above  0°K.  This  has  been  recently  shown  by  Betts  and 
( 2 ) 

Elliott  j  see  also  Section  3  in  this  thesis.  Therefore 

it  would  be  very  interesting  indeed,  to  know  the  exact 

nature  of  the  transition  point  and  any  singularities  which 

occur  there  and  to  compare  these  with  other  three  dimensional 

loose  packed  antif erromagnet  lattices.  So  far  any  available 

information  has  come  from  high  temperature  series  expansions, 

but  all  these  attempts  do  not  produce  any  conclusive  facts 

( 3 ) 

about  low  temperature  behaviour.  For  example  Burley  , 

using  the  Kikuchi  approximat ion;  extrapolated  series 

expansions  for  various  lattices  including  the  f.c.c.  into 

( 4 ) 

the  low  temperature  region.  Sykes  and  Fisher  ,  using 

improved  techniques,  have  extrapolated  exact  series 

expansions  of  the  loose  packed  sq  and  HC . (honeycomb )  anti- 

(  S  6  )  / 

ferromagnet  lattices.  BakerW3  '  found  the  Pade  approximant 

method  very  successful  in  extrapolating  ant if erromagnet 

susceptibility  series  expansions  of  almost  all  the  loose 

( 7 ) 

packed  and  close  packed  lattices.  In  1961  Danielian 
showed  that  two  dimensional  order  existed  in  the  f.c.c. 
antiferromagnet  lattice  and  he  calculated  the  ground  state 
entropy.  Later  in  1964  ^ 8 ^  he  calculated  the  first  three  terms 
of  the  low  temperature  partition  function  expansion  and  gave 
a  rough  estimate  of  the  transition  temperature. 
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5. 


1.3  Scope  of  Thesis 

Section  2  outlines  the  development  of  the 
exact  series  expansions  which  have  been  applied  In  sub¬ 
sequent  sections.  C.  Domb^  and  Newell  and  Montroll ^ 10 ^ 
have  contributed  excellent  review  articles  which  cover 
most  of  the  material  in  this  section. 

In  section  3, the  low  temperature  series  expansions 
of  the  partition  function  and  a  long-range  order  parameter  are 
derived  for  the  f.c.c.  ant if erromagnet .  A  discussion  of  the 
configurational  problems  involved  is  also  included. 

In  section  4,  the  critical  temperature  is  estimated 
by  extrapolating  series  expansions  of  the  specific  heat 
curve  and  the  long  range  order  parameter.  The  critical 
entropy  is  also  estimated.  Section  5  is  a  brief  summary 
of  the  results  obtained  and  conclusions. 
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Theory 


2.1  Low  Temperature  Series  Expansions 


The  partition  function  can  be  expressed 


N 


N1’N12 


g(N5NiaN12)  exp(-E(N,N1JN12)/kT) 


(3) 


where  the  conblnatorial  factor  g(N3NlsN12)  specifies 
the  number  of  possible  configurations  of  the  lattice 
corresponding  to  atoms  with  spin  up  and  N^2  nearest 
neighbour  bonds  between  unlike  spins.  E  is  the  total 
energy  of  the  microstate  and  depends  only  on  N,N^  and  N^2 

An  equivalent  low  temperature  expansion  of  (3) 
can  be  formed  by  starting  with  the  lattice  in  the  ground 
state  and  enumerating  all  possible  configurations  as  more 
and  more  spins  have  overturned. 

Thus  we  can  write 


ZN  =  H  g(N.a.b)  exp  (-E(N,a,b)/kT)  (4) 

a ,  b 

where  a  =  number  of  spins  overturned 

bJ  =  increase  in  energy  resulting  from  the  a  overturned 
spins . 

g(NjO,o)  is  then  the  degeneracy  of  the  ground  state  con¬ 
figuration  . 


7. 


We  first  consider  the  ferromagnetic  case 
(J>o).  The  total  configurational  energy  E  can  be  written 


E  =  E  +  E, 
o  b 


where  E  the  ground  state  energy  is 

N 


E  =  -N  m  H  -  q  — —  J 
o  2 


q  =  the  number  of  nearest  neighbours  per  atom. 


E,  =  bJ  +ea  m  H 
b 


Substituting  in  (4)  we  get 


z>r=z-q(N/4)  p-(N/2)  j-  g(N,a)b)  zb/2  f/2 


N 


b 


(5) 


We  have  introduced  the  low  temperature  variables 


z  =  exp (-2  J/kT) 


y  =  exp— ( 2mH/kT ) 


which  tend  to  zero  as  T  0. 

If  we  write  the  series  in  (5)  as  TV  N(y  ,z  )  then  (  5  )  Is 


shortened  to 


'N 


=  z-Q(n/4)  h-(n/3)  yyN(U)Z) 


(6) 


In  zero  field  y  =  1  and 


-q (N/4 )  W  (l  z) 

z  tr ±>Z) 


1 


8. 


For  the  loose  packed  antiferromagnetic  lattices 
in  zero  field,  there  is  an  exact  correspondence  in  energy 
between  the  configurations  of  a  ferromagnet  and  those  of 
the  ant  if erromagnet .  This  means  that  the  partition  function 
in  the  absence  of  a  magnetic  field  must  be  invariant  under 
the  transformation  J  -*  -J  (see  for  example  (9)  page  175)* 
However,  we  will  primarily  be  concerned  with 
the  close  packed  f.c.c.  antiferromagnetic  lattice  which 
does  not  fall  into  this  category.  From  here  on 
all  statements  will  refer  to  the  f.c.c.  Ising  antif erromagnet 
unless  otherwise  indicated.  Its  partition  function  is 
not  invariant  under  such  a  transformation. 

Let  us  suppose  that  in  the  ground  state  in  zero 
field  every  +1  spin  is  along  the  z  direction  (direction  1) 
and  every  -1  spin  is  antiparallel  to  the  z  direction  (direction 
2).  We  will  further  suppose  that  every  +1  spin  will  have 
q  nearest  neighbours  which  are  also  +1  spin  and  q-^ 
nearest  neighbours  which  are  spin  -1  (the  subscripts 
referring  to  the  direction  of  each  spin). 

The  following  relations  are  then  true 

qll  +  q12  =  q22  +  q21  =  q 
q12  =  q21 
q22  =  qll 


* 
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9. 


It  also  follows  that  there  will  be  N/2  +1 
spins  and  N/2  -1  spins.  The  above  conditions  are  not 
necessarily  true  for  all  close  packed  lattices.  The  two 
dimensional  triangular  lattice  is  an  example  in  which 
they  are  not  ^ ^ . 

In  a  strong  magnetic  field  the  ground  state  must 

become  f erromagnet ically  ordered  but  for  weak  fields  the 

ground  state  will  be  unchanged. 

The  critical  field  H  occurs  when  these  two 

c 

states  have  equal  energy  i.e. 


qN 


(-J)  -  mN  H  = 
v  '  c 


N 


2 


(-Jqn  +  q12  J  ) 


(8) 


Prom  equation  (8)  we  can  determine  the  critical  field 


mH  =  -q, 0  J 
c  12 


(J<o)  (9) 

We  will  now  consider  the  low  temperature  partition  function 


expansion  for  fields  less  than  H( 
We  have 


E  = 
o 


NJ 


2 

=  NJ  q  ’ 


'(q12  qll) 


where  q ’ 


=  a  -  -£L 

q12  2 


=  -bJ  -  2m  H  (ab  -  a2) 


. 


10. 


where 


=  no.  of  1  spins  overturned 
=  no.  of  2  spins  overturned 


N 


2:  2:  g(N,a,b)  exp[(-Nq’  J  +  bJ  +  2m  H(a. -ap ) )/kT] 

a,a^  b  d 


Z 


N 


(10) 


a,  a1  b 


where  y  =  exp(2J/kT)  is  the  antiferromagnetic  low  temperature 
variable  which  tends  to  zero  as  T  -+  0 .  We  will  make  use  of  (10) 
later.  As  an  illustration,  the  term  corresponding  to  two  over¬ 
turned  spins;  one  in  direction  1  and  the  other  in  direction  2 
is 


(11) 


2.2  High  Temperature  Series  Expansions* 


Starting  from  equation  (2)  in  the  form 

z„=  H  ft"  exp (k  c . a . )  /iT  exp  ( 

N  o1=fl  nn  1  J  i 


it  can  be  shown  that  this  reduces  to 


£=o  m=o 


(12) 


*  see  references  (9)  and  (10) 
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where  A  =  number  of  nearest  neighbour  bonds 
K  =  J/kT 
u  =  tanh  K 
I  =  mH/kT 
v  =  tanh  I 

G(£,m,N)  is  the  combinatorial  factor  corres¬ 
ponding  to  the  number  of  graphs  of  £  bonds  and  m  distinct 
spins;  the  m  spins  coinciding  with  the  odd  vertices.  An 
odd  vertex  is  a  vertex  at  which  an  odd  number  of  bonds 


meet . 


FIGURE  I 


t 


■-© - & 


m  =  2 

£  =  5 
(a) 


m  =  2 

£  =  13 

(b ) 


m  =  1 

£  =  6 

(c) 


In  Figure  l,in  which  the  spins  are  circled,  (a) 
and  (b)  contribute  but  (c)  does  not. 

The  expansion  (12)  is  a  high  temperature  series 
expansion,  since  as  T  ->  00  u  and  v  -*  0 .  (12)  is  valid  for 
an  ant  if  erromagnet .  Under  the  transformation  J  -»■  -J 


[u(J)] 


(-if  [u(+J)f 


r.rf  i  **  0  ^  >c  *  i  '  o  j  1  'H oq 


i 


c  -h 
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In  zero  field  the  partition  function  reduces 
to  a  series  involving  closed  polygons  only. 

2.3  Thermodynamic  Properties 

The  thermodynamic  properties  are  obtained 
directly  from  the  series  expansions  of  the  partition 
function.  Some  of  the  thermodynamic  variables  are: 

The  internal  energy  per  spin 

U  =  kT2  — - ln  (13) 

9T  N 

The  specific  heat 

Cv  =  (9U/9T)y 

=  — (kT2  — )  (ln  Z  )/N  (14) 

9  T  9  T 

The  entropy  per  spin 

rT  C 

S  =J  dT  ~y~  (15) 

The  magnetic  properties  can  also  be  obtained.  The 
magnetization  per  spin  along  the  z  direction  is  the  average 
induced  magnetic  moment  per  spin  and  is  given  by 

9  In  Z]\f 

9  H  N 


M  =  kT 


(16) 


T3i  •»  U 


i  ;s..  e  ■  2  ef  T 


TP 
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The  susceptibility  per  spin  Is 

X  =  3M/3H  (17) 

It  is  more  convenient  to  express  these  variables  in  terms 
of  y  and  y .  Using  the  relation 


T 


3 

3  T 


-(log  y)y  — 


the  thermodynamic  variables  become 


U  =  -  2J  y-fy  (m  Zn)/N 

a 

-~-=  (log  y)2  y  [  3  /  3  y  +  y  ( 3  2/3y2 )  ]  ( log  ZN/N) 
S  \  [Cv/(y  In  y)  ]  dy 


M  =  -2m  y(3/3y)  (log  ZN/N) 


(2m) 

X  =  - — - 

kT 


y[3/3y  +  y ( 3 2/3 y 2 ) ] ( log  ZN/N) 


The  zero  field  susceptibility. 


lim 
H^O  X 


(18) 

(19) 

(20) 
(21) 
(22) 
(23) 


(24) 


;  s . 


; 


3.  Applications  to  the  P.C.C.  Antif erromagnet 


3.1  Introduction 

An  exact  low  temperature  series  expansion  for 
the  partition  function  and  a  low  temperature  expansion  for 
a  long  range  order  parameter  has  been  calculated  by  con¬ 
figuration  counting  methods. 

( 7 ) 

Danielian  has  completely  classified  the 

ground  state,  showing  that  it  was  highly  degenerate  but 

the  ground  state  entropy  per  spin  was  zero.  He  showed  that 

two  dimensional  ordering  in  one  of  three  perpendicular 

(2) 

planes  must  exist.  Recently  it  was  shown  that 

the  two  dimensional  order  can  be  expressed  by  a  three 
dimensional  long  range  order  parameter  averaged  over  the 
statistical  ensemble.  The  transition  temperature  has  been 
estimated  by  finding  the  temperature  at  which  the  order 
disappears . 

3.2  The  Ground  State 

(7 ) 

According  to  Danielian  the  nature  of  the  ground 

state  can  be  elucidated  as  follows: 

The  f.c.c.  lattice  of  N  sites  can  be  subdivided 
into  N  tetrahedra  such  that  each  bond  is  contained  in  one 
and  only  one  tetrahedron. 


asnrlq 
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Thus  there  are  a  total  of  6N  bonds.  In  the 
ground  state  each  tetrahedron  must  be  in  Its  own  ground 
state;  I.e.  two  spins  must  be  up  and  two  down.  Thus  the 
configurational  ground  state  energy  is  +  2NJ.  The 
degeneracy  of  the  ground  state  is  seen  most  easily  by 
considering  the  lattice  as  being  built  up  of  layers 
of  two  dimensional  quadratic  lattices  so  that  each  atom 
of  one  layer  is  located  above  the  middle  of  a  square 
in  the  lower  layer  and  each  atom  in  the  lower  square 
is  a  nearest  neighbour  of  the  atom  centered  above  that  square. 
This  is  illustrated  in  Figure  2. 

F.C.C.  ANTIFERROMAGNETIC  GROUND  STATE 


(110)  direction 

- - - - ^ 

lower  layer  bond 

- upper 

-  bonds  between 

adjacent  layers 


«  -  spin  up 
o  -  spin  down 


FIGURE  2 
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A  partial  layer  of  25  atoms  is  shown  and  9 
atoms  of  the  next  layer.  If  we  begin  construction  of  the 
lattice  with  the  lower  layer  antif erromagnet ically  ordered 
(alternating  black  and  white)  then  when  we  arbitrarily  add 
an  atom  to  the  next  layer  its  spin  state  completely  deter¬ 
mines  the  states  of  all  the  atoms  of  that  layer.  Referring 
to  Figure  2,  if  we  choose  atom  1  of  the  upper  layer  black  then 
atom  2  must  be  white  since  atoms  1,2,3  and  4  form  a  tetra¬ 
hedron  in  which  we  can  only  have  2  black  atoms. 

We  conclude  that  the  ground  state  consists  of  layers 
of  antif erromagnetically  ordered  quadratic  nets,  also  each 
atom  has  4  like  nearest  neighbours  and  8  unlike  nearest 

I/O 

neighbours.  There  are  (2N)  such  layers  and  so  the 


(2N) 


1/3 


degeneracy  of  the  ground  state  is  in  the  order  of  2 
since  there  are  two  possible  states  for  each  ordered  plane 


The  ground  state  entropy  is  then 


S ( o )  =  k(2N)1/3  In  2 

The  ground  state  entropy  per  spin  tends  to  zero 
in  the  limit  N  +  «. 


3.3.  Long  Range  Order 

We  designate  two  sites  a,b  along  a  line  connecting  nearest 
neighbour  atoms  in  the  lattice  and  specify  that  the  distance  between 
them  be  an  odd  number  of  bond  lengths.  If  the  expectation  <  a  > 
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Is  non  zero  we  say  that  order  exists  and  If  <  a  a,  >  is 

a  b 

non  zero  in  the  limit  distance  a  -  b  -*  00  then  long  range 
order  is  present. 

For  the  f.c.c.  ant if erromagnet  we  know  that  in 

the  ensemble  of  states  at  T  =  0°K  one  third  of  all  possible 

a  and  b  sites  will  be  in  the  ordered  planes.  For  these 

cases  a  a,  =  -1.  The  expectation  <  a  a,  >  for  all  other 
a  d  a  d 

cases  we  would  expect  to  be  zero  from  the  way  in  which  the 
ground  state  is  constructed.  The  expectation  <  c  a,  >  for 
any  a,  b  then  is  not  zero  but  exactly  -1/3.  We  can 
conclude  then,  that  "partial"  three  dimensional  order  exists. 

The  parameter  a  is  not  symmetrical  in  three 
dimensions  since  we  do  not  know  if  a,  and  b  lie  in  the 
ordered  plane  or  in  one  of  the  other  two.  However,  we 
can  construct  an  order  parameter  that  is  symmetrical 
in  all  3  directions.  Consider  a  large  regular  tetrahedron 
in  the  lattice  built  up  of  smaller  tetrahedra,  and  with 
its  edges  parallel  to  the  nearest  neighbour  directions. 

If  each  edge  contains  an  even  number  of  spins,  the  following 
property  is  evident.  In  the  ground  state  a  pair  of  opposite 
edges  of  the  tetrahedron  will  lie  on  parallel  planes. 

Since  there  are  three  pairs  of  opposite  edges,  corresponding 
to  the  three  sets  of  mutually  perpendicular  planes,  one  of 
the  pairs  must  lie  in  the  ordered  plane.  As  a  consequence 


18 . 


of  the  alternation  of  spins  along  these  edges,  two  of 
vertices  will  have  "spin  up"  and  two  "spin  down". 

Let  <  >  represent  the  mean  fraction  of 

tetrahedra  of  a  given  size  with  two  spins  up  and  two 
spins  down.  Also  define 

<  ai  >  =  mean  fraction  of  spins  up 

<  a2  >  =  mean  fraction  of  spins  down. 

Clearly,  then  if  there  is  no  correlation  among  the  spins 
of  the  4  vertices  of  the  tetrahedron  we  must  have 

2  2 

<t22  >  ^  6  <  >  <  a  2  > 


(24) 


since  there  are  6  distinct  ways  of  labelling  the  4  spins. 
Following  the  method  of  Domb^^  we  modify  (24)  to  take  care 
of  the  correlation  by  introducing  an  order  parameter  R^  as 
follows 

2  2  5 

c  1 2  2  >  =  6  <  a  i  >  <  a2  >  ( 1  +  ^  )  (25) 

If  =  0  there  is  no  correlation.  For  a  zero 
field  configuration  <  aj  >  =  <  a2  >  =  1/2  so  we  can  write 


5/3  (Rjj )  =  8/3  (<  t22  >)  -  1  (26) 

and  we  see  that  long  range  order  disappears  when  <  t22  > 

decreases  to  3/8. 
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In  the  following  equationsR^  is  expressed  in 

terms  of  the  Ising  variables  a  .a,  3a  ,  and  a  where  a3b3c 

a  d  c  Q 

and  d  denote  the  4  vertices  of  the  tetrahedron. 


t 
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F  "  F  (oa°b  +  °b° 

+  rf  (a  a,  a  a  , ) 
o  a  b  c  d' 


+  a  a  ,  + 

c  d 


a  a  +  a  a  ,  + 

a  c  ad 


a  a 
b 


(27  a) 


Also 


R4  ' 


-  ^  (a  a,  +  a,  a  +  a  a +  a  a  +  c  a,  +  cna  ,  ) 
5  a  b  be  cd  ac  ad  bd 


+  5  (oaabac°d) 


(27  b) 


The  right  side  of  (27  a)  is  always  unity  for  two  plus  and  two 
minus  spins  and  zero  otherwise.  Equation  (26)  will  be  used 
later  to  determine  the  Curie  point  defined  as  the  temperature 
at  which  ->  0 . 

It  was  mentioned  earlier  that  is  equal  to 

-  1/3  in  the  ground  state  and  equal  to  zero  at  the  transition 
temperature.  This  suggests  a  natural  definition  of  a  two 
spin  order  parameter. 

R2=-3<a1a2>  (28) 

Although  R2  lacks  the  symmetric  features  of  R^ 
it  can  still  be  thought  of  as  a  three  dimensional  order 
parameter  as  long  as  we  do  not  specify  the  planes  in  which 


. 
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sites  1  and  2  lie.  In  fact  it  is  possible  to  show  that 
and  R^  are  related  by  the  equation 

5R4  =  2R2  +  3(R2)2  (29) 


Equation  (29)  is  derived  from  a  relation  between  the  Ising 

variables  <  a a ^  >  for  2  spins  and  <  >  for  4  spins. 

We  have  2 

<  >  =  1/9  <  °^a2  since  two  of  the  four 

spins  of  the  tetrahedron  are  independent  of  the  spin  states 

of  the  other  two. 

We  can  see  from  (29)  that  in  the  range  0  £ 
we  have  R2  the  equality  being  true  only  at  the  end 

points.  Physically  this  means  that  the  correlation  between 
two  spins  is  always  greater  than  the  correlation  between  four 
spins  which  is  what  one  would  expect. 

For  the  calculation  of  Tc ,  the  parameter  R^  was 
used  although  R2  would  lead  to  the  same  result. 


■ 
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3.4  Low  Temperature  Series  Expansions 


The  first  few  terms  of  the  expansion  series 


equation  (10) 

are  listed  in 

and  their  corresponding  graphs 

Table  I,  below: 

TABLE  I 

and  contributions 

TERM 

GRAPH 

CONTRIBUTION 

g  ( N  3  o ,  o ) 

1 

g (N, 1, 8 ) 

•»o 

N  ,  -Is 

2  (vi  +  t  ) 

g (N , 2 , 12 ) 

• - 0 

4  N 

g(Ns2,l6) 

0  0  J 

p  p  p 

•  • 

(1/4)(N  -10N)(y  +U 

•  0 

(1/4)(N2  -  16N ) 

g(N, 3 , 16 ) 

o— m— o 

ION ( y +y-1 ) 

g ( N , 4 , 16 ) 

• — 9 

O - 0 

C  3/2 )N 

LOWEST  ENERGY  OVERTURNED  SPIN  CONFIGURATIONS 


. 
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The  configuration  increase  in  energy  -bJ  is  given  by 

AE  =  -8nJ  -  4J  6  *  (30) 

n 

where  n  =  number  of  spins  in  the  configuration 

6  =  number  of  even  bonds  less  the  number  of  odd  bonds 

n 

(an  even  bond  connects  like  spins  and  an  odd  bond  connects 
unlike  spins ) . 

The  smallest  change  in  energy  is  obtained  by 

overturning  one  spin  and  is  equal  to  -8J.  This  gives  rise 

4  - 1  4 

to  the  factor  (N/2)u  y  for  an  overturned  black  spin  and  (N/2)y  y 

for  an  overturned  white  spin.  The  next  smallest  change  of 

energy  possible  is  -6Ja  produced  by  overturning  adjacent 

black  and  white  spins.  The  number  of  configurations  increases 

rapidly  for  high  powers  of  y.  All  the  ant  if erromagnet  low 

temperature  configurations  for  the  f.c.c.  up  to  5  spins  and 

all  others  corresponding  to  an  energy  change  of  kjE  =  14J  or 

less  have  been  enumerated  in  Appendix  A.  Table  II  on  the 

next  page  shows  how  the  number  of  topologically  different 

graphs  increases  with  AE . 

*  (30)  follows  directly  from  the  ground  state  conditions  that 
every  spin  has  4  like  neighbours  and  8  unlike  neighbours . 
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TABLE  II 


AE/2 

Number  of  Distinct  Graphs 

0  J 

0 

4  J 

1 

6  J 

1 

8  J 

3 

10  J 

8 

12  J 

31 

14  J 

127 

NUMBER  OF  OVERTURNED  SPIN  CONFIGURATIONS 
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3.5  Configuration  Counting 

The  number  of  ways  a  graph  can  be  placed  on  a 
lattice  of  N  sites  divided  by  N  is  called  the  "lattice 
constant"  for  that  graph.  The  enumeration  of  all  the  graphs 

and  the  calculation  of  the  lattice  constants  is  a  difficult 

(  q  ) 

and  time  consuming  task,  Domb,  in  his  review  article  has 
outlined  some  of  the  standard  techniques  for  counting,  there¬ 
fore,  only  the  difficulties  encountered  with  antiferromagnetic 
f.c.c.  configurations  will  be  pointed  out  here. 

The  first  problem  in  obtaining  the  low  temperature 
series  expansion  is  to  enumerate  all  the  graphs  corresponding 
to  each  increment  of  energy.  This  is  easy  enough  for  low 
energy  graphs  but  is  liable  to  error  for  higher  energy  graphs 
as  a  glance  at  Table  II. will  indicate.  One  approach  is  to 
assume  that  all  the  graphs  have  been  accounted  for  when  one 
can  not  think  of  any  more.  However,  equation  (30)  seems  to 
offer  a  more  systematic  approach. 

For  example,  suppose  we  want  to  enumerate  all 
possible  6  spin  configurations  that  produce  an  energy  change 
of  AE  =  28J  (excluding  any  magnetic  contribution)  when  the 
spins  are  overturned.  Equation  (30)  states 

«6  =  +  7  - 


2  =  5 


?jBri 

- 
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i.e.  The  number  of  odd  bonds  must  be  five  greater  than 
the  number  of  even  bonds.  So  the  only  possibilities  are 
(5a0),  (6,1),  (7,2),  ....  etc.  where  the  numbers  in 
brackets  are  the  numbers  of  odd  and  even  bonds  respectively. 
We  begin  with  (5,0)  and  we  notice  that  the  only  possible 
colourings  for  the  f.c.c.  are  2b,  4w;  2w,  4b;  and  3b,  3w. 

The  first  two  colourings  are  equivalent.  We  consider  each 
colouring  separately  and  group  all  of  our  graphs  accordingly. 

In  the  group  G[(5,0)  3b,  3w]  we  find  the  following 
set  of  graphs: 


(1) 

o— •— o 

•  c>— 4 


(2) 

o  4—o 


(3) 


(4) 


(5) 


o— »■ 


o 


-o- 


1 


-o- 


#— o- 


(6) 


(7) 


O— ®  a? 

4—0  6 


Figure  3 


and  we  wish  to  know  if  it  is  complete. 

If  the  6  vertices  are  all  considered  to  be 
distinguishable  (numbered  for  instance),  then  the  number  of 
ways  of  adding  the  5  odd  bonds  is  clearly  or  126. 


. 
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Now  define  the  order  P  of  each  graph  as 


P  = 


n,  !  n 
b  w 


where  =  number  of  black  vertices 

nr  =  number  of  white  vertices 
w 

S  is  a  symmetry  number  which  is  equal  to  the  number  of 
automorphisms  of  the  labeled  graph. 

In  our  example 


 (3 


=  I! 

2 

 (3 


 (3 


4 

(3 

T 


)2 


1 


i2 

)2 


) 


=  18 
=  18 
=  36 

=  9 
=  9 


i 


The  total  2j_  =  126  showing  that  all  the  graphs  in  group 
G[(5j0)  3b  3w]  have  been  accounted  for. 

The  only  discouraging  feature  of  this  method  is 
that  some  configurations  that  are  incompatible  with  the  lattice 
structure  must  also  be  included. 

The  task  of  determining  the  lattice  constants  for 
each  configuration  is  complicated  by  the  degeneracy  of  the 
ground  state.  For  example,  the  number  of  times  the  con¬ 
figuration  P  T  appears  in  the  lattice  is  dependent  upon 


which  particular  ground  state  of  the  ensemble  of  ground  states 
that  the  lattice  is  in.  Consequently  it  is  necessary  to  use 


elementary  probability  considerations. 


w  %. 

■ 

« 
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3.6  Long  Range  Order  Series  Expansion 

Returning  to  equation  (26),  which  defines  the 
long  range  order  parameter  R^,  we  can  compute  the  expectation 
<  t^2  >  by  taking  an  ensemble  average  as  follows: 

■5T  g(N,b  jtpp )  yb/2 

<  t22  >  =  .  .  (3D 

g  ( N , b  )  yb/2 

where 

if  t22  =  0 
if  t22  =  1 

Thus  g(N3bat22)  counts  the  number  of  ways  each  configuration 

can  be  placed  on  the  lattice  so  that  the  left  side  of  (27) 

is  one  after  overturning  the  spins.  Clearly  the  only  possibilities 

are 

1)  none  of  the  4  vertices  are  overturned 

2)  one  w  and  one  b  vertex  are  overturned 

3)  two  w  and  two  b  vertices  are  overturned 

The  four  vertices  are  assumed  to  be  infinitely  far  apart  so 
that  no  closed  configurations  can  overlap  two  vertices 
simultaneously . 

3.7  Counting  Theorem 

A  relationship  between  the  terms  in  equation  (31) 
to  the  terms  in  the  partition  function  has  been  rigorously 


g(N,b  a t22 )  = 


0 


g(Nab) 


•  .. 
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established,  and  considerably  simplifies  computation. 

We  will  adopt  the  following  notation: 

Let  X^  represent  a  connected  configuration 

{X^}  -  the  contribution  of  X  to  the  coefficient 
in  the  partition  function 
[X.]  -  the  contribution  of  X  to  the  coefficient 
in  the  expansion  series  for  (8/5)  R^. 
XfXjX^.  -  a  configuration  with  three  components. 


{XXX2  —  X.  —  xM} 
CX1X2  Xi  XM 


(x1> 


[X.] 


M 

M 


are  similar  contributions  from  configurations  of  M 
components . 
o 

X^  -  a  graph  which  coincides  with  one  of  the  two  white  vertices 
of  the  tetrahedron.  We  will  refer  to  such  graphs  as  being 
"rooted'1 . 

Thus 

o  ®  @ 

[X.  X.  X,  is  the  number  of  ways  the  disconnected 
i  j  k  M  J 

graph  of  M  components  can  be  put  on  the  lattice  so  that  two 

black  vertices  and  one  white  vertex  of  the  tetrahedron  are 

overturned.  Note  that  it  is  impossible  for  one  connected  graph 

to  be  rooted  to  two  vertices  simultaneously  since  the  vertices 

are  assumed  to  be  infinitely  far  apart. 


* 


GCt:  O  Li 


* 


; 

* 
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The  counting  theorem  then  states  for  a  disconnected  configuration 
of  M  distinct  components 


[Xi]M  =  {x/}m  “  (VN)  n,  {X,} 


'i  M 


i  L  iJ  M 


o  {X. }  {X.} 

“p  {Xi>M  - 1 - ^  (w  W  +B.B.  +  2W  B.  )X 

N  ij  {X.  X.  }  1  J  1  J  1  J  ij 

1  J 


3  {Xi}M 


(X  }{X  }{X  } 

- - - *2— — —  (W.  W.B,  +B.  W.W,  )  X 


+ 


N 
384 


ij5k  {x.xjxk} 


i  j  k  i  j  k  i  j  k 


N 


W 


{X4  } 


i  M  / 


(X.  }{X.  HX,  }{X.  } 

— - - *2 - - - -  ( W.  W  .  B,  B  )X  .  ..  n 

1  j  k  '  ijk£ 


{xixjW 


(32) 


The  corresponding  counting  theorem  for  the  coefficients  of  the 
expansion  series  of  (2/3)R2  is 


[X.]M  =  {Xi}M  -  (2/N)  n.  {X.}M 


o  {X . } {X  .  > 

+  (4/tr)  {x.  }M  E  -J - (W.  W.+B.B.  )  X  .  . 

1  i,J  {X.X.}  1  J  1  J  1J 


(33) 


1  J 


where  W.  =  the  number  of  white  vertices  of  X. 

1  1 

B.  =  the  number  of  black  vertices  of  X. 

J  ..  J 

n.  =  the  number  of  vertices  for  X. 

1  1 

The  summations  in  (32)  and  (33)  are  taken  over 
all  different  combinations  of  the  M  components  so  that  in  each 
term  i  ^  j  ^  k  ^  &  .  The  X’s  are  factors  to  take  care  of. cases 
in  which  the  M  components  are  not  all  distinct.  The  proof  of 
the  counting  theorem  is  given  in  Appendix  B  along  with  the 
values  of  the  X  factors. 
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4.  Thermodynamic  Properties 


4.1  Low  Temperature  Series  Expansion  of  Thermodynamic 

Variables 


Applying  the  methods  described  in  sections  2 
and  3j  the  following  low  temperature  expansion  for  the 
partition  function  in  a  magnetic  field  was  obtained. 

i/3  ' 

ZN  =2(2N)  y-N  {1  +  (N/2)  (y+y_1)y4  +  4Ny6  +  [(N2/4)  -  5/2  N 
+  ION  (y+y-2)  +  ((N2/8)  -  ( 5/4 ) )  (y2+y-2 )  ]y8 

p  _ ip  _p  i n 

+[ (209/2 )N  +  (2N  -27N)(y+y  )  +  9N(y  +y  ) ]y 

+[18N2  -  (619/2)  +  (1/16)(N3  -  30N2  +  6735N) (y+y-1) 

+(5N2  -  85N)(y2+y-2)  +  (1/48)(N3  -  30N2  +  320N ) ( y 3+y "3 ) ]y 12 
+(1/8 ) [8N3  -  296N2  +  38265N  +  (774  N2  -  16049  N) (y+y-1) 

+  (4N3  -  148N2  +  7290N)(y2  +  y-2)  +  (36N2  -  768N) ( y 3  +  y "3 )  ]y 14 


+ 


} 


(34) 
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On  physical  grounds 

zn  ^  Cz(y)]N  where  Z(y)  is  the 
partition  function  per  spin, 

log  Z  =  N^«  ni/N)  log  ZN(y,N)]  (35: 

We  get 

log  Z  =  -log  y  +  { (1/2)  (y+y-1)  y4  +  4y6 

+[-(5/2)  +  10  (y+y  Z)  -  (5/4)*  (y2+y  2)]y3 

+[(209/2)  -  27  (y+y-1)  +  9(y2+y"2)]y10 

+  [-(619/2)  +  (6735/16 )  (y+y-1)  -  85(y2+y-2) 

+  (20/3)  (y3+y-3)]y12  +  ( 1/8 ) [38265  -  16049  (y+y-1) 

+  7290  (y2+y'2)  -  768  (y 3+v"3 ) ]y14  +  - }  (36) 


*  Equation  5  in  reference  8  is  in  error  for  this  coefficient. 


.  . 


<  ;  "  '•  ■  '  } 


32 . 


In  zero  field  the  partition  function  reduces  to 

Z  =  -log  y  +  y4  +  4y8  +  15y8  +  (137/2)  y10 

+(9017/24)  y12  +  (19211/8)  y14  +  -  (37) 

From  (36)  and  (37)  we  derive  the  low  temperature  series 
expansions  for  the  thermodynamic  variables  as  described  in 
Section  2. 

Internal  Energy  per  Spin 

U  =  2J{  1  -  4y4  -  24y6  -  120y8  -  685  y10  -  (901^/2)  y12 

-  (134477/4)  y14 - }  (38) 

Zero  Field  Specific  Heat 

(C  /k)  =  (log  y)2  (iSy11  +  l44y^  +  960y^  +  6850y10  + 

+  54l02y12  +  (941339/2)  y14  + - )  (39) 

M  =  -2m  y-N  ( (N/2 ) ( u-y"1 ) y 4 

+ [ ION ( y-  v1 )  +  ((N2/4)  -  (5/2)N)(u2-u"2)]y8 


Magnetization 
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+[(2N2  -  27N)(u-ti_1)  +  l8N(u2-u"2)]y10 
+[(l/l6)(N3  -  30N2+6735N) (p-u-1)+(10N2-170N) (u2-y 
+ ( 1/1 6 ) (N3  -  30N2  +  320N)(y3-y~3)]y12 
+(l/8)[(774N2-l6049N)(y-y-1)+(  8N3-296N2+l4580N ) 
(y2-y_2)+(108N2  -  2304N)(y3-y_3)]yll4+ - } 

Magnetization  per  Spin 

M  =  -2m  {(l/2)(y-y_1)y4 

+[10(y-y_1)  -  ( 5/2 )  (y 2  -y  '2)]y8 

+[-27 (y-y-1)  +  l8(y2  -y"2)]y10 

+[(6735/16) (y-y-1)  -  170(y2-y”2)+20(y3-y'3)]y12 

+  ( 1/8)  [-16049  (y-y_1)  +  l4580(y2-y'2)-230it(y3-y"3)]y 14 


(40) 


} 


(41) 


t  •  1  l  ; 


■ 


i  ■  •  4,' 01  * 
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Susceptibility  per  Spin 

x  =  ^  [(2m)2/kT]  {  ( 1/ 2  )  ( y+y""1 )  y4 

+[10(y+y  1)  -  5 (y2+y  2)]y^ 

+[-27(y+y~1)  +  36  (y2fy~2)]yl0 

+[ (6735/16 ) (y+y  2)  -340(y‘~  +  y  2 )  +  60  (y2+y  2)]y12 
+(l/8)[-l6049(y+y  2)  +  29160 (y^  +  y  2 ) -6912 (y -  +y  ^)]y24 
+ - }  (42) 

taking  the  limit  - 

X  0  =  [(2m)2/kT]  {y4+10y8+l8y10+(2255/8)  y  12+  ( 6199/4  )y14+  — }( 43  ) 

The  two  spin  long  range  order  parameter  series  expansion  is 
R2  =  1  -  (l/3)[12y4+96y6+624y8+4l88y10+(60585/2)y12 

+  235059  y14 - ]  (44) 

The  4  spin  long  range  order  parameter  is  given  by 

R4  =  1  -  (8/5)(4y4+32y6+202y8+1300y10+(l8l79/2)y12+6946lyl4+— ) 

(45) 


■ 


{ - 4 


» 
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4.2  Series  Extrapolations  Method 

To  locate  the  singularities  in  the  functions 

whose  series  have  been  listed  in  the  previous  section, 

( q ) 

the  ratio  method  of  Domb  and  Sykes v  and  the  more  recent 
Pade  approximant  method  have  been  used.  Both  methods  apply 
when  the  terms  of  the  series  in  question  have  reached  some 
asymtotic  behaviour.  In  references  (5) and  (6),  Baker  has  given 
a  detailed  account  of  the  applicability  of  the  Fade  approximant 
method,  A  brief  explanation  of  the  two  methods  follows  but 
for  further  details  the  above  mentioned  references  should  be 
consulted , 


4*2-1  Ratio  Method 


Assuming  an  infinite  power  series  has  the  asymtotic 
form  of  a  function  (1  -  x/x  )Y  near  some  point  xp  on  the  x 

W  V 


axis 


i .  e . 


00 

(1-  */x  )Y  ^  21  axn. 

n=sD 


I,  . 

Then  the  limiting  form  of  the-  ratio  of  the  coefficients  is 
given  by 


n+1 

a 


n 


*  ““  [1  -  (  Y/n+l) 

c 


( 1/n+l ) ] 


(46) 


. 

-  • 
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Therefore,  if  the  ratio  (an+1)/(an)  is  plotted  against 

l/(n+l)  the  relation  should  become  linear  for  sufficiently 

large  n.  The  intercept  on  the  (a  +^)/(an)  axis  will  be 

1/x  and  the  limiting  slope  will  be  -(y+l)/(x  ). 
c  c 

Figures  4  and  5  show  a  plot  of  the  ratios  of  the 
series  for  the  long  range  order  parameter  and  the  ratios 

of  the  specific  heat  series  respectively.  From  previous 

* 

studies  the  long  range  order  parameter  is  expected  to  behave 

asymtotically  like  [1-(T/T  )]2e  near  T  or  in  terms  of  y, 

like  [l-(y2/y  2)]2^*  0  is  (1/8)  for  the  two  dimensional 

0 

lattices  and  is  approximately  (5/16)  for  three  dimensional 
lattices.  By  extrapolating  the  curve  of  the  ratios  in  Figure 
4  our  estimates  are 

kT 

— -  -  1.72  1.12 
J 

6  =  .4  ±.l 

The  specific  heat  curve  singularity  at  T  is  possibly  of 
the  form 

\  (1  -  (T/T  ))"1/b  (47) 

k  c 

^  (l  -  (y  /yc  ) ) 

as  was  indicated  for  the  ferromagnetic  f.c.c.  and  triangular 

( Q ) 

lattices  by  Domb  and  Sykes'7  .  It  is  also  possible  that  the 


* 


. 


. 
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singularity  is  logarithmically  infinite >  corresponding  to 

( 1/b )  =  0,  on  the  low  temperature  side  of  T  ,  characteristic 

0 

of  2  dimensional  lattices.  If  the  singularity  is  logarithmically 
infinite  then  the  asymtotic  form  of  the  ratios  of  the  coefficients 
of  the  series  [ (d/dy) (Cv/k) ] should  be  constant. 

The  ratios  are  plotted  in  Figure  5(b)  and  show  a 
slight  increase  for  large  n.-  Again  there  does  not  seem  to  be 
enough  terms  to  give  conclusive  evidence  of  the  asymtotic  be¬ 
haviour  of  the  series. 

An  extrapolation  of  the  ratios  of  the  series  of 
equation  (47)  failed  to  give  a  positive  value  for  1/b 
(Figure  5(b))  but  the  critical  point  was  estimated  to  be 

kTc 


J 


1.6  ±  .1 
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.  < 

4.2.2  Pade  Approximant  Method 

If  a  function  is  known  to  possess  a  singularity 

of  the  type  f(x)  =  (x  -  x  then  its  logarithmic 

0 

derivative  will  Joave  a  simple  pole  and  residue  equal  to  -y  at  xc 

i.e.  — —  [In  f(x)]  =  "■  Y| 

dx  x-x 

c 

In  the  Pade  approximant  method  the  series  representing  the 
logarithmic  derivative  of  the  function  is  equated  to  a  ratio 
of  two  polynomials  PM(x)/Q^(x).  The  numerator  is  of  degree 
M  and  the  denominator  of  degree  N3  where  we  choose  N  +  M  +1 
equal  to  the  highest  power  in  our  power  series.  The  unknown 
coefficients  of  P^Cx)  and  Q^(x)  are  determined  by  equating 
coefficients  of  like  powers  of  x  and  solving  N  +  M  linear 
equations.  The  ratio  P^(x)/Q^(x)  is  called  the  [N,  M]  Pad§ 
approximant  and  we  look  for  the  pole  of  (d/dx)[ln  f(x)] 
among  the  N  roots  of  Q^(x)  and  y  will  be  found  by  calculating 
the  residue  at  the  pole. 

Long  Range  Order  Approximants 

Table  III  lists  values  of  Tc  and  3  from  Pade 

2 

approximant  analysis  of  ( d/du ) ln[R2 ]  where  u=  y  .  The 
approximants  have  been  listed  in  Appendix  C. 


» 
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TABLE  III 


( kTc  /  J ) 

Pade  Approximant 

3 

[i,i] 

1.85 

.072 

[2,2] 

1.81 

.064 

[1,3] 

1.85 

.069 

[3,1] 

- 

«** 

[2,33 

1.93 

.092 

1.57 

.0016 

[3,2] 

1.725 

,021 

CM] 

1.84 

.049 

[4,1] 

- 

[5,0] 

- 

- 

FADE  APPROXIMANT  ANALYI 

» 

IIS  OF  LONG  RANGE  ORDER  2 

3ERIES' 

#  No  real  roots  found  ntar  Tc . 


i 


g(S  ' 


4o . 


Specific  Heat 

In  Table  IV,  estimates  of  Tc  and  1/b  were  obtained 

_p 

by  analysing  the  logarithmic  derivative  of  y  (log  y)  Cv/k. 

4  2 

The  factor  y  (log  y)  was  removed  from  the  series  expansion 
since  it  does  not  appreciably  influence  the  shape  of  the 
curve  near  the  sigularity. 


TABLE  IV 


Pade  Approximant 

( kTe  /J ) 

1/b 

[4,0] 

no  real  roots 

near  Tc 

[3,1] 

1.96 

.246 

[2,2] 

1.65 

.235 

[1,3] 

1.66 

.168 

PADE  APPROXIMANT  RESULTS  OF  THE  SPECIFIC  HEAT  SINGULARITY 


V  0  r-  s  •  ,  1  3,-  '  '^9 


•  ■ 

• 
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Magnetic  Susceptibility 

It  is  well  known  that  ferromagnetic  susceptibilities 
are  infinite  at  the  transition  temperature  but  for  anti¬ 
ferromagnetic  lattices  the  behaviour  is  quite  different.  For 

the  loose  packed  quadratic  and  honeycomb  lattices,  Sykes  and 
(in 

Fisherv  have  confirmed  that  the  susceptibility  is  finite, 
exhibiting  a  vertical  slope  at  the  transition  point  and  then 
reaching  a  maximum  slightly  above  Tc .  Baker ^ ^  and  Burley 
have  analysed  the  high  temperature  antiferromagnetic  sus¬ 
ceptibility  series  but  their  results  become  uncertain  at  low 
temperatures . 

It  is  reasonable  to  assume  then  that  the  magnetic 
susceptibility  of  the  f.c.c.  antiferromagnetic  lattice  behaves 
like  the  loose  packed  antiferromagnetic  lattices  near  the 
transition  point.  If  so,  we  should  expect  a  sigularity  (pole) 
in  the  derivative  of  xq  and  possibly  a  "zero”  above  Tc .  Table 
V  lists  the  results  of  a  Pade*  approximant  analysis  of  the 
derivative  of  the  susceptibility  series.  The  series  so  far 
seems  quite  irregular  and  ratios  of  the  coefficients  do  not 
indicate  any  asymtotic  behaviour.  However,  a  Pade*  approximant 
analysis  of  the  derivative  of  the  susceptibility  series 
was  attempted  using  two  approximant s .  The  results  of  the 


. 

•  , 

'  . 

>i 

« 
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calculations  are  listed  below  in  Table  V.  The  [1,4]  approxi- 
mant  gave  a  somewhat  higher  answer  for  kTc/J  than  previous 
estimates  and  the  [2,3]  approximant  failed  completely. 


TABLE  V 


/ 

Pade  Approximant 

kTc/J 

[1,4] 

2.15 

[2,3] 

- 

/ 

PADE  APPROXIMANT  ANALYSIS  OF  THE 

SUSCEPTIBILITY  SERIES 

* 

• 
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4.3  Critical  Entropy 

The  entropy  can  be  calculated  using  equation  (15) 
and  directly  integrating  the  specific  heat  curve.  It  is 
of  interest  to  know  the  critical  entropy  given  by 


For  all  lattices  is  k  In  2  =  0.6931  k.  The  percentage 

increase  of  the  entropy  to  T  increases  with  coordination 

c 

number.  For  the  f.c.c.  ferromagnetic  lattice  (q  =  12) 

S  /k  =  0.597  so  that  a  very  large  portion  of  the  entropy 

change  takes  place  below  T  .  For  the  two  dimensional 

c 

quadratic  lattice  S  /k  =  0.306  is  much  less.  Our  estimate 
for  the  f.c.c.  antiferromagnetic  lattice  is  Sc/k  =  0.312 
which  is  remarkably  near  to  the  value  for  the  quadratic 
lattice.  This  is  strong  evidence  that  the  effective 
dimensionality  of  the  f.c.c.  antiferromagnetic  lattice  is 
two  rather  than  three. 

The  integral  in  equation  (48)  was  obtained  by  fitting 

the  specific  heat  series  to  a  logarithmic  singularity  at 

( 12 ) 

T  ,  similar  to  a  method  used  by  Wakefield 
c 

We  equate  the  specific  heat  series  to  a  logarithmically 


infinite  function  as  follows. 


vqoior  '  ;  l  inn 


H  3C3  v<  J  3bT  ‘f  t  •  f  ^  ;  r-; 


hi 


.  • 1  r' 3  r.  ri :  •  ;  c  t  r 


£  D^n.^i  ;b30  j,  cl  *  oi«xr  ••?  --a  o.-Vroeq-i  & 


•»''•>  ■'■  ■:  •  »r  '  ■  ,  :  .ti  i f ; .. 
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=  [log  y]2  £l6y4  +  l^^y6  +  960y8  +  6850y10  +  54l02y12 

i  ii  *) 

+  (9i»l339/2)y  + - ( 

O 

7 

=  A  ( In  y)2  In  (1  -  (y2/yc2))  +  (In  y)2  0-2  y2n  (49) 

n=o 


The  first  term  on  the  right  side  of  (49)  represents  the 
singularity  and  the  polynomial  ensures  that  function  well 
represents  the  specific  heat  for  values  of  y  much  less  than 
yc .  We  chose  yc  =  .1307  from  the  [1,3]  Pade  Approximant 
in  Table  IV.  Equating  coefficients  of  equal  powers  of  y 
we  end  up  with  eight  equations  and  nine  unknowns,  A,  a  , 

- ai4*  However,  we  can  set  a^  equal  to  zero  and 

solve  for  the  remaining  eight  unknown  coefficients.  Setting 
a^  to  zero  is  effectively  saying  that  our  series  has  settled 
down  to  asymtotic  behaviour.  The  new  series  obtained  in 
this  way  is: 


=  2.13802  (In  y)2  In  [1  -  ( y 2/ ( . 13069 ) 2  )  ] 

-  (In  y)2  (16.35932  y2  +  46.58828  y4 

+  175.2723  y6  +  872.238  y8  +  4365.785  y10 
+  17414.454  y12}  (50) 


V 


,0 


‘ 
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In  terms  of  the  variable  x  =  T/Tc ,  the  critical 
entropy  is  given  by 


k 


In  Figure  6  (Cv/kx)  is  plotted  against  x  using  equation 


(50).  Thus  the  critical  entropy  is  given  by  the  area 
under  the  curve  and  we  find  Sc/k  =  .312  for  the  f.c.c. 


We  have  also  plotted  the  specific  heat  curves 


for  the  antiferromagnetic  quadratic  lattice  and  the 
ferromagnetic  f.c.c.  lattice  for  comparison.  The  area 
under  the  curves  for  the  f.c.c.  ant  if erromagnet  and  the 
simple  quadratic  lattice  are  almost  identical  (.312  and  .306 
respectively)  but  their  shapes  differ  significantly.  We 
see  also  that  a  much  greater  percentage  of  the  entropy  change 
occurs  below  Tc  for  the  f.c.c.  ferromagnet  than  for  the 
corresponding  antiferromagnet . 


■  ■ 


APPLICATION  OF  RATIO  METHOD  TO  LONG  RANGE  ORDER  SERIES 


The  ratios 

7 

where 


a=0 


(a2n+2^/^a2n^  are  Pitted  against  l/n 

a2n  y2n  "  =  R2  =  1  ‘  4y4  “  32y6  '  208y8 

-  I396y10  -  ( 20195/2 )y12  -  78353yl4 
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APPLICATION  OF  RATIO  METHOD  TO  SPECIFIC  HEAT 

SERIES 


•j— i  •  c/  \  l/(n+l) 

Figure  5(a) 


The  ratios  (a2n-t2^/^a2n^  are  P10^6^  against  l/n 

7 

where  21  a~  y2n  =  l6y  +  l44y  +  9^0y  +  6880y10 

n=2 

+  54l02y12  +( 941339/2 )y14 


APPLICATION  OP  RATIO  METHOD  TO  THE  DERIVATIVE  OF 


THE  SPECIFIC  HEAT  SERIES 


The  ratios  (a 


6 

z 

n=l 


w 


2n 


2n+2//a2n^  are  Plot,ted  against  l/(n+l)  wher 

d  [  i _  £v  -j 

d(y2)  (log  y)2  k 

32y2  +  i*32y4  +  3480y6  +  34,250y8 
3246l2y10  +  3294686 . 5y12 


8v?s„  +  yQS^Z  +  y££U  +  s^££ 


COMPARISON  OF  SPECIFIC  HEAT  CURVES  FOR  SIMPLE 
QUADRATIC  (EXACT),  F.C.C.  FERROMAGNETIC^  AND 
F . C , C .  ANTIFERROMAGNETIC . 
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5.  Summary  and  Discussion 


5 . 1  Summary 


The  problem  of  determining  the  low  temperature 
properties  of  the  f.c.c.  Ising  model  antiferromagnetic 
lattice  has  long  been  an  outstanding  problem  in  theoretical 
physics,  mostly  because  of  the  complicated  ground  state 
configuration.  The  only  feasible  approach  to  the  problem 
is  by  configuration  counting  methods  which  quickly  be¬ 
comes  an  enormous  task  as  more  and  more  terms  in  the 
expansion  series  are  required.  Previously  only  three 


terms  had  been  calculated 


(8) 


but  we  have  now  extended  the 


series  to  seven  terms  and  many  of  the  contributions  to 
higher  terms  have  also  been  calculated  (see  Appendix  A). 

An  extensive  study  of  long  range  order  in  the  f.c.c. 
lattice  was  also  made.  A  series  representation  for  a 
three  dimensional  order  parameter  which  measures  the 
correlation  among  four  spins  at  the  vertices  of  a  large 
tetrahedron  was  calculated  and  a  relation  between  this 
four  spin  correlation  and  the  more  usual  two  spin  correlation 
parameter  was  obtained.  To  obtain  the  long  range  order  ex¬ 
pansion  series  special  counting  theorems  were  developed 
to  simplify  the  calculations. 


. 
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The  transition  temperature  Tc  is  characterized 

by  the  disappearance  of  long  range  order.  Using  the  Pad 

(56)  ( 9 ) 

approximant  method  5  and  the  ratio  method  we 

have  estimated  that  the  long  range  order  disappears  at 

— =  1.9  ±  .1 

J 

From  the  low  temperature  expansion  series  of  the 
partition  function  in  a  magnetic  field  H  series  of  many 
other  thermodynamic  functions  have  been  calculated 
including  the  internal  energy  U,  the  specific  heat  C  , 
the  entropy  S,  the  magnetization  M  and  the  zero  field 
susceptibility  xo* 

Since  the  specific  heat  series  was  somewhat 
regular  in  comparison  to  the  series  for  the  other  thermo¬ 
dynamic  functions,  and  also  because  it  is  expected  to 
exhibit  a  weak  singularity  at  the  transition  point,  our 
studies  were  concentrated  on  its  properties. 

In  general  the  Pade  approximant  and  ratio  method 
of  series  extrapolation  confirmed  the  existence  of  a 
transition  point  near  (kTc/J)  =  1.9  and  although  there 
is  still  doubt,  the  singularity  in  Cv  seems  to  be 
logarithmic  in  nature  on  the  low  temperature  side  of  . 


CDS 


. 
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The  specific  heat  series  was  fitted  to  a  logarithmic 
singularity  below  T  and  compared  to  the  specific  heat 
curves  of  the  two  dimensional  quadratic  and  the  f.c.c. 
ferromagnetic  lattices  (Figure  6).  It  was  interesting 
that  the  curve  for  the  f.c.c.  antif erromagnet  more  closely 
resembled  the  curve  for  the  quadratic  lattice  than  that  of  the 
f.c.c.  ferromagnet.  Possible  reasons  for  this  are  discussed 
in  the  next  section. 

The  critical  entropy  S  was  calculated  by  direct 
integration  of  the  specific  heat  curve. 

Tc 

i.e.  sc  =  J  (Cv/T)dT 

The  estimate  for  S  was  0.312k  as  compared  to  0.306k  for 
the  simple  quadratic  lattice  and  0.597k  for  the  f.c.c.  anti¬ 
ferromagnetic  lattice,.  The  implications  here  are  dis¬ 
cussed  in  the  next  section. 

The  susceptibility  series  was  too  irregular  to 
attempt  an  extrapolation.  It  appears  at  least  two  more 
terms  will  be  required  before  the  behaviour  near  T  can 
be  determined. 

For  low  temperature  in  the  range  T/T  <  .5  all 

w 

the  series  expansions  appear  to  converge  rapidly  and 
therefore  the  thermodynamic  functions  are  well  represented. 
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by  the  first  few  terms  of  their  associated  series.  But 
in  general  it  was  difficult  to  determine  the  asymtotic 

i 

form  of  the  series  because  of  this  insufficient  number  of 
terms.  However,  the  graphs  of  the  coefficient  ratios 
of  the  specific  heat  and  the  long  range  order  series 

indicate  that  the  coefficients  may  be  approaching  asymtotic 
behaviour  and  if  this  is  the  case?  two  or  three  more  terms 

t 

should  be  sufficient  to  get  some  conclusive  answers. 
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5.2  Effective  Dimensionality. 

The  ground  state  has  been  shown  to  consist  of  layers 
of  antif erromagnetically  ordered  quadratic  nets  in  which 
the  spin  states  do  not  seem  to  influence  the  spin  states  in 
adjacent  ordered  layers,,  at  least  for  very  low  temperatures. 
Because  of  this  unusual  ordering  in  the  ground  state  the  question 
naturally  arose  as  to  whether  the  thermodynamic  properties 
would  be  characteristic  of  the  two  dimensional  quadratic 
antiferromagnetic  lattice  or  would  they  more  closely  resemble 
the  properties  of  the  three  dimensional  close  packed  ferromagnetic 
lattices . 

The  answer  is  probably  something  in  between.  From 
the  point  of  view  of  order  among  the  spins  in  the  lattice 
sites  we  have  shown  that  order  is  essentially  only  in  two 
dimens ions. But  we  can  still  define  three  dimensional  order 
parameters  which  correlate  two  and  four  spins  at  arbitrary 
locations  on  the  lattice  as  long  as  they  lie  along  a  straight 
line  of  odd  number  of  bond  lengths  for  two  spins  and  are  at  the 
vertices  of  a  tetrahedron  for  four  spins.  The  specific  heat 
curves  (figure  6)  show  that  there  is  a  remarkable  difference 
between  the  specific  heat  of  f.c.c.  antif erromagnet  and  the 
f.c.c.  f erromagnet;  a  much  larger  portion  of  the  curve  for 
f.c.c.  f erromagnet  being  below  Tc.  However  the  f.c.c.  anti¬ 
ferromagnetic  curve  is  similar  to  that  of  the  simple  quadratic 


■ 


. 

. 
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except  for  the  difference  in  area  distribution.  This  difference 
must  somehow  be  related  to  the  way  in  which  the  quadratic 
layers  of  the  f.c.c.  antiferromagnetic  lattice  interact, 
otherwise  the  curves  would  be  identical  in  nature. 

Our  estimate  is  Sc  =  ,312k  showing  that  about  45^  of 
the  total  entopy  change  takes  place  below  Tc .  For  the  ferro¬ 
magnetic  and  loose  packed  lattices  Sc  is  roughly  proportional 
to  q.  Sc  =.306k  for  the  simple  quadratic  (q  =  4)  and  Sc=.597k 
for  the  f.c.c.  ferromagnetic  lattice  (q  =  12).  However 
our  estimate  for  the  f.c.c.  antif erromagnet,  being  much 
lower  than  the  critical  entropy  for  the  f.c.c.  ferromagnet 
seems  to  indicate  that  it  is  not  really  q  that  decides  the 
magnitude  of  Sc  but  possible  an  "effective  q"  where  in  general 
the  "effective  q"  would  be  less  than  q  for  close  packed  anti¬ 
ferromagnetic  lattices  and  equal  to  q  for  ferromagnetic  and 
antiferromagnetic  loose  packed  lattices. 


. 
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5.3  Comments  on  Experimental  Measurements 

On  past  occasions,  the  Ising  model  has  produced 
reasonable  good  qualitative  agreement  with  experiment 
and  we  should  expect  the  f.c.c.  antif erromagnet  lattice 
to  be  no  exception.  There  are  many  substances  in 
nature  which  have  a  f.c.c.  antiferromagnetic  structure. 

Among  them  are  MnO,  MnS,  FeO,  CoO,  NiO,  and  the  salts 
K^IrClg  and  (NH^ ) 2lrClg .  MnO  has  been  studied  by 
Bizette*  et  al  and  exhibits  a  logarithmic  type  singularity 
in  the  specific  heat.  The  susceptibility  reaches  a 

maximum  at  the  critical  point. 

( 13 ) 

Domb  and  Miedema  have  written  an  excellent 
summary  of  experimental  measurements  of  transitions  in 
solids  which  include  a  table**  of  thermodynamic  properties 
of  antiferromagnetic  salts  with  loose  packed  structures. 

For  the  cases  of  spin  1/2,  the  values  of  S  /k  are  smaller 
than  the  corresponding  close  packed  antiferromagnetic  values 
which  seem  to  agree  with  our  result . 


*  See  reference  (15),  p  436  for  graphs  of  MnO  properties 


**  See  Table  9,  p  340  of  reference  (13) 


r'  f  edtcf Iris  bn*  e  *  *9*;f9S.t0 
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5.4  Further  Studies. 

There  are  some  properties  of  the  f.c.c.  spin  1/2 
Ising  model  which  have  not  been  studied  in  this  thesis  but 
are  still  of  considerable  interest  to  theoretical  physicists. 
There  is  the  calculation  of  middle  range  order  parameters 
(between  second,  third  etc.,  nearest  neighbours).  The  short 
range  order,  between  nearest  neighbours  is  linearly  related 
to  the  internal  energy  and  therefore  has  essentially  been 
calculated.  A  study  of  all  these  parameters  would  more  clearly 
specify  the  manner  in  which  order  spreads  as  temperature  is 
decreased  to  0°. 

The  critical  energy  would  also  be  of  interest  since 
this  quantity  can  be  compared  to  experimental  measurements. 

Very  little  can  be  said  about  the  susceptibility 
series  at  present  because  of  irregularity  in  the  first  few 
terms  of  the  series.  Therefore  more  terms  will  be  required 
before  information  about  the  behaviour  near  the  transition 
temperature  and  especially  about  the  maximum  above  Tc  can 
be  obtained. 


. 
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APPENDIX  A 


Face  Centered  Cubic  Antiferromagnetic  Low 

Temperature  Lattice  Configurations  * 


^Notation:  The  superscript  of  "C"  is  half  the  energy  associated 

with  overturning  the  spins.  The  subscript  numbers 

the  configurations  and.  specifies  the  number  of 

22 

vertices  in  the  graph.  Thus,  Is  associated 

with  energy  44J,  has  five  vertices  and  is  the 
96th  five  spin  graph. 
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~  (10  N*-  57  8  N  4-  90  24) 
o 

24 

C  0-0-0  0  o 

r/6/ 

(6N2-414  N  +  7  7  72) 

24 

C  Oo©  0  o 

ZUL 

~(3Na-  2  2  0  N  4-  4248) 

C 0  0-0-0  Q  o 

i76J 

-g  (16  N2- 1 1 1  6  N  +21  512) 

c;*  o-o-o  o  o 

.17  b* 

±(20  N -1  3  8  6  N  4-  26308) 

8  ' 

£  o-o-o  o  © 

J-(  8  N2-—  5  2  8  N  -f-  9264) 

8 

16 

C  0-00  o  © 
>v*6 

±(10  N*—  688  N  +12648) 

8 

20 

C  O  ©  0  ©  © 

J~/6? 

C°  0  0  0  G  0 

r/6t j 

C  °  o  o  O  O  0 

■ 
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13 
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n 

S*l 


AE  /2 

O 


=  12J 
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H 

c 

0 - 0 


/4 

c 


Q — Q — Q — O 


6 — 6 — o 


75 


/4 

c 

g// 


11 

4 


c 

SlS 


■O 


6—6 — 6 


21 


/  4 


O — O — O — O 


819  0—0 — 0 — 0 


48 
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2 


/4 

C 

9*9 


C 
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1_9 

4 


7 


1_3 
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1Z 
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2 


35 
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8 


12  Vertices  AE/2  =  l4J 


11 

8 
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APPENDIX  B. 


Counting  Theorem 


To  prove  equation  (32)  we  proceed  as  follows: 
Consider  a  graph  X^  X^  --  X^  --  X^  of  M  components 
and  the  number  of  ways  that  this  graph  can  be  placed  on 
the  lattice  so  that  two  vertices  of  the  tetrahedron  are 
black  and  two  are  white  is 


[X„]M  =  «.}„  - 


(X,  } 


i  M 


o  o 

(X,X0 } 
1  J 


M 


•  • 

txiV 


M 


o 


-  {xixjVm 


•  O  O 

(X.X.X,  }  M 
l  j  k  M 


o  o  «  • 

-  {xixjW 


M 


( B  .  1 ) 


o 

Now  first  consider  {X^}^ 

The  number  of  ways  X^  can  be  rooted  to  a  white 
vertex  of  the  tetrahedron  Is 


2  {X.}  W. 

i  l 

{o} 


—  U,  }  W 
N 


The  factor  2  occurs  because  there  are  2  white  vertices  on 
the  tetrahedron;  because  each  of  the  white  vertices 


' 


i"  ?n  •  ■ 


. 

. 
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{x±} 

of  X.  could  be  the  root  and  there  are  — - —  ways  of 
1  (o> 

orienting  X,  on  the  root.  However,  we  are  not  rooting 

{X^}  but  {X1>M.  The  number  of  ways  this  can  be  done  is 


4 


N 


{Xj_}  wi 


<X1>M 

{Xi} 


W 


lXi}M 


But  in  doing  this  some  other  of  the  M  components  are 
going  to  be  rooted  to  the  other  three  vertices.  We  must 
subtract  these.  So  we  get 


o  • 


o  o 


{ Xi }  M 


N 


Wi  {xi}m  -  s:  <xiV 


M 


-  2  21  {X±X  } 

ij  J 


~  2 


o  o 


o  #  tf 


'0  0  •  0 


i  J 


{XiXJXk>M 


{XiXJXk>M  -  2 


i  ,  j  a  k ,  £ 


<xixjWM 

(B.2) 


Similarly 


{Xi  jM 


4 


•  o 


•  • 


N 


B 


i  {xfM  -  £  {xiVm  -  2  £  {xixj} 
bJ  /-j 


M 


-  2X  (X,*  X  }  -  ZlfX.X  X  }  .  2  Z  a±x  XkXt) 

|,J,K  J  ;j,K  J  0 


M 


(B.3) 
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Substitute  B.2  and  B03  into  B.l 


ExiN  -  (x1]M 


7T~  {Vm  (Bi+Wi)  +  2  ^  <*iV« 

IM  '  j  J 


O  O  O  «  0  0# 

+  2Zltxixj}M  +2  21{X1X1}m  +  3'^{X1XiXk} 


»,J 


i  j  M 


fjj 


,jJ;K 


i  j  k  M 


+  3.Z  +  te/hXjW  M 


(B .  4 ) 


•  » 


Now  consider  {X^Xj}^ 


Clearly 


8  « 


{xixj> 


{x1){xj} 

i  i  2 
;  8  1 


2  B^Bj 


8 


N 


«  {X.HX.}  B .  B. 

2  i  J  i  J 


•  «  o  { X »  }  m 

.7  {X  X  }  =  -42-  (X  }{X  }  B  B  — 

l  J  n  N  l  j  l  J  {xix,} 

1  J 


®  •  O  0 


ZL  {X,  x,  x  xi  . 

•uZl  '  J  K  xrA 


(B.5) 


^  Xi  }]yj 

since  - — -  is  the  number  of  ways  of  putting  the  M-2 

{W 


configurations  on  the  lattice  after  fixing  X.  and  X 


-  •  '  -  ■■  : 


' 
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similarly 


o  o 

(x  X  ^ 

1  j  M 


8 


<X.  > 


i  M 


*X.  Hx.  >  W.  W. 

N  i  J  i  J  {XoX} 

i  J 


—.00®  00®® 

-  2.  tx.x.x  -5”  {X.X.X,X„>M 

...  i  j  k  M  <L —  l  j  k  z  M 

i,  j,  Kj  t 


(B.6) 


and 


U.X  >  =  -W-  (x.HxJ  W.B 


VM 

i“j  M  2"  lAiJ‘V  wi°j  {  } 

1  J 


_  o  •  o  0  0® 

2.I-{XXX}  -2^(XXX) 

ilJ.K  j  VI;  K  j 


_ _  O  •  O  9 

2.  (x,x,x  xt)  M 


(B  .  7  ) 


Now  substituting  equations  (B.5)*  (B.6)  and  (B.7)  into  (B.4) 


Cx.]M  -  {Xl}M 


4 


N 


{X±}M  (B.  +  W.) 


o  ^  {X.}{X.} 

+  { X .  >  — ±—2—  (W.W.  +  B.B.  +  2  W.B.) 

l  M  r„  „  ,  i  J  i  J  l  J 


N 


2 


{X.  X . } 


»jJ  '"1  "j 


\ 


O'  o  • 


3Z_  «lXjXk}M 

WJ;  K 


-  3  2I{xi*jVM  -6.S 

i  j  J;  K  ij  Jj  K 


0  0®# 

{X  X  X  X  } 


M 


(B.8) 


■ 


. 
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Now  continuing  in  the  same  fashion. 


oo® 


IjOj  K 


{xiVk}M 


2x2x8 


N 


3 


W.WoB,  - 

i  3  k 


{X1HX  }{XkHXi}]y[ 

{xihxk} 


(B .  9 ) 


-  2 


o  o 
{X.X 


1 1  £ 


IVI 


0  •  t 

*>.  tX.X.X.  } 
~  i  j  k 

»  ji 


32 


N 


-  W.W.W, 

3  i  j  k 


{XjLHXJHXkHX1}M 

txiW 


oo® 


2  „£,(xiTOjm 


UJj  KjJL 


Substituting  (B.9)  and  (B.10)  into  (B.8) 


N 


{Xi}M  <Bi  +  Wl> 


{X.HX.} 
{X,  },-  >  . - . 


n2  l"1JM  br  «jxj} 


(WjWj  +  +  2W1Bj) 
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{X.  HX,  }{X,  > 


«.}M  !E— — L— —  (w,w 

N3  «1XJXk} 


i'M,rr  7„  +  h'f'V 


-5T 


Q  ©  ®  • 


+  6,2 

ij  J,  i<,  4 
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The  last  term  in  (B.ll)  is 


__  o  o  •  • 

>  {X.X.X,  Xn  }  M 

i  j  k  £  M 

j  ->j  Kl 


2x2xl6 


{xi}{xiHxkHX2 

{xixj W 


{ Xj_  Im 


x  W .  W  o  B,  B 
l  j  k 


(B . 12 ) 


Finally  we  arrive  at  the  result 


[  x  i  ]  M 


N 


ni  ^i^M 


{X  }{X  } 

+  — {X,h„  Z_  — (W.W.  +  B.B.  +  2W.B.)  X. 


N 


1  M  £j~  {X.X.}  '  1  J  _1  J  1  1J 

1  J 


96 


N 


3 


{X,  }{X,HX.  } 

i  J  W 


(W.  W.B,  +B.  W.W,  )  X  . 


'i  M  f-—  (XXXI  1  J  k  i  j  k'  ijk 
04 ^  Ai  j  k 


+ 


384 


{Xi}M 


{xpu .  ){XkHX£} 


iw  {xixixkx2 


( W .  W  .  B,  B  )  X  .  . ,  /  T3 

i  J  k  ijk  (B. 


Where  we  have  replaced  B^  +  W^  by 


13) 


. 
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The  X  factors  have  been  included  to  take  care  of  identical 
components  in  the  M  component  graph.  It  can  be  shown 


that 


X.  i 

1 


Xijk  1 


Xi  *  Xj  *  Xk 


=  3 _  x 

Ai 


Xj  Xl  i  Xk 


V 


Xijk£  1 

5 

"  “£■ 
1 

“  “T 

l 

TE 


x.  i  Xj  f  xk  i  x 

(2  of  X's  identical) 
(3  XTs  identical) 

(4  X’s  identical) 


For  M  =  1  only  the  first  two  terms  in  (B.13)  are  non  zero 
and  the  calculation  involved  is  trivial. 

For  M  =  2  the  third  term  becomes  non  zero  but  it  can  be 


simplified  as  follows: 
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[X±X2 ]  =  (X1X2) 


4n 


N 


{x1x2> 


8 


N 


2  {X1HX2>  [n1n2  +  3(W1B2  +  B^)]  X12 


The  corresponding  counting  theorem  for  the  two  spin 
correlation  parameter  R2  is  of -course,  simpler  in  form 


and  is 


CXl]M  “  tXl}M  "  N  ni  fXl}M 


n  U,HX,} 

+  — * —  { X  .  }  M  . . . —  *■—  (W.B.  +  Bi  W.  ) 

1  11  i ,  J  {XiXJ  }  1  J  1  J 


(B.15) 
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APPENDIX  C. 


Fade  Approximants 


C.l  Long  Range  Order  Approximants 


d  log  (R2) 
d(y2) 


8y2  -  96y 4  -  8 N Oy 6  -  7l40y8 


62585y10  -  57l459y12 


[1,1]  = 


-  8y  2  -  26y4 

l  -  8 . 75y2 


[2  2]  -  -  8y 2  -  100.1905/  -  258 . 6i9yb 

1  +  . 5238095y2  -~78 . 95833yir 


[1,31 


8y2  -  25.87672ylj  +  1.47936y6  +  222.9444y8 

1  -  8 .  7654ly2 


[3,1] 


-  8y  ^  +  ,l6912yi4 


1  -  12 . 021l4y  -  39.25368y  +  840.26386y 


F 


[2  3]  -  ~  8y 2  +  70.697625y4  +  339.4893y6  +  512.6643y8 

1  -  20.8372y2  +  102.6l027yir 


[3,2] 


-  8y 2  +  54.8391y^  +  l40.8803y6 

1  -  18 .85^8 9y ^  +  103.64863/”-  I56.52025y^" 
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8y  ^  +  6717  _ 

1  -  851  •  70588yr'  + . lOii^ThlSSy-  32489. 0294y6  +  84388 . 3l6l8y8 

8y2  -  22.95259y4  +  36.568890y6  +  529-977790y8 
_  +  2609 . 8ll2l67y1Q  _ 

L  l  -  9.l309259y2 


[5,0] 


12y2  +  3900y4 


8y2  _ 

100. 5y6  -  2.i25y8  -  1784. 375y10 


C  .  2 


Specific  Heat  Approximants 


d 


In 


dy 


1 


Cv 


(log  y)2  l6y  k 


=  9  +  120y2  +  393 . 375y4  +  3792y6  +  ^42224. 47y: 


[4,0]  = 


1  -  13.3333y2  +  139-06944y4  -  I626.l48l5y6  +  16748 . 19334y: 


[3,1] 


_ 9  +  92.69373y2  _ 

1  -  3 . 0340299y2  -  3 . 25460295y!|  -  245 . 3262475y 6 
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ro  9  +  12.8480136yd  -  968 . 460579y 4 

L  ^ ^  J  p  ~TT~~ 

1  -  11.9057766y  +  7 - 42862^0y 

n  ni  _  9  +  19* 78369^ 6y 2  -  942,8424y^  -  588.28768y6 
Uj  jj  p 

1  -  Ii.i35l45y 


C.3 


Zero  Field  Susceptibility  Approximants 


kT 


dy  \ ( 2mr ) 


X0=  2(2y3  +  40y^  +  90y9  +  l691.25yl:L 


+  10848. 25y13  + 


[1,4]  =  y3 


2-12 ,828677y2+40 . POOPy11-^ . 5735403y6+ni3 ■ 95953y8 

1  -  6.4l43385y2 


[2,3]  =  y3 


2  -  9-46l5449y2  -  23.274024y4  -  99.230898y6 

1  -  4.7307725y2  -  31 . 63701196y3 
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